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ABSTRACT 

We report that neutral hydrogen (H i) gas clouds, resembling High Velocity Clouds (HVCs) observed 
in the Milky Way (MW), appear in MW-sized disk galaxies formed in high-resolution Lambda Cold 
Dark Matter (ACDM) cosmological simulations which include gas-dynamics, radiative cooling, star 
formation, supernova feedback, and metal enrichment. Two such disk galaxies are analyzed, and Hi 
column density and velocity distributions in all-sky Aitoff projections are constructed. The simulations 
demonstrate that ACDM is able to create galaxies with sufficient numbers of anomalous velocity 
gas clouds consistent with the HVCs observed within the MW, and that they are found within a 
galactocentric radius of 150 kpc. We also find that one of the galaxies has a polar gas ring, with radius 
30 kpc, which appears as a large structure of HVCs in the Aitoff projection. Such large structures 
may share an origin similar to extended HVCs observed in the MW, such as Complex C. 
Subject headings: methods: N-body simulations - galaxies: formation - galaxies: ISM - Galaxy: 
evolution - Galaxy: halo - ISM: kinematics and dynamics 



1. INTRODUCTION 

High Velocity Clouds (HVCs) are neutral hydrogen 
(Hi) gas clo uds with velocities inconsistent with galac- 
tic rotation ( Wakke ~ fc van Woerdenlll997ft . From our 
vantage point within the Galaxy, they appear to cover a 
large portion of the sky relatively isotropically. HVCs 
do not appea r to possess a stellar component (e.g. 
IPimon Blitz| UMl IPiegel (£jp HI) and their dis- 
tances (and masses) are generally unknown. Direct dis- 
tance constrai nts have only been made for a sel ect num- 
ber of HVCs (|WakkeiH200ll: iThom et al.ll200l . There 
are still open questions as to whether HVCs are local 
to the Milky Way (MW) or distributed throughout the 
Local Group (LG); whether they are peculiar to the 
MW or are common in disk galaxies; whether they are 
gravitationally bound or pressure confined; whether they 
contain Dark Matter (DM); and their degree of metal- 
e nrichment. 

iPisano et alJ (|2004T) report that there are no HVC-like 
objects with H I mass in excess of 4 x 10 5 M in three LG- 
analogs. They suggest that if HVCs are a generic feature, 
they must be clustere d within 160 kpc of the host galaxy, 
ruling out the original blitz et all (|l999f) model in which 
HVCs are gas cloud s distributed in filaments on Mpc- 
scales. Thilk er et all (|2004l) and IWestmeier et al.1 ([2005) 
find 16 HVCs around M31, with Hi masses ranging from 
10 4 to 6 x 10 5 M . Most of the HVCs are at a projected 
distance < 15 kpc from the disk of M31. Some clouds ap- 
pear to be gravitationally dominated by either DM or as 
yet undetected ionized gas. They also found two popula- 
tions of clouds, with some of the HVCs appearing to be 
part of a tidal stream, and others appearing to be primor- 
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dial DM dominated clouds, left over from the formation 
of the LG. Finally, it has also been speculated that cool- 
ing instabilities in Cold Dark Matter (CDM) halos could 
lead to clouds within 150 kpc of MW type ga laxies, with 
HVC-like properties ( Mai ler fc Bullocldl2004fl . 

In this Letter, we report that these mysterious Hi 
clouds also appear in MW-size disk galaxies formed in 
ACDM cosmological simulations. We demonstrate that 
the simulated galaxies show HVCs comparable in pop- 
ulation to the observed ones. We also find that large 
HVCs resembling Complex C appear in simulated galax- 
ies. Therefore, we conclude that HVCs appear to be 
a natural byproduct of galaxy formation in the ACDM 
Universe. The next section describes our methodology, 
including a brief description of the simulations and how 
we "observe" the simulated disk galaxies. In Sectional 
we show our results and discuss our findings. 

2. METHODOLOGY 

We analyze two disk galaxy models found in cos- 
mological simulations that use the multi-mass tech- 
nique to self-consistently model the large-scale tidal 
field, while simulating the galactic disk at high resolu- 
tion. These simulations include self-consistently many 
of the the important physical processes in galaxy forma- 
tion, such as self-gravity, hydrodynamics, radiative cool- 
ing, star formation, supernova feedback, and metal en- 
richment. The disk galaxi es we analyz e corre spond to 
"KGCD" and "AGCD" in iBailin et, al . (2005), and we 
use these names hereafter. Both simulations are carried 
out with our galactic ch emodynamics package GCD+ 
((Kawata fc Gibsonll2003h . 

The deta ils of these simulations are given in 

Baili n et alJ ((20051) . Table H summarizes the simulation 
parameters and the properties of the galaxies. Column 1 
is the galaxy name; Column 2, the virial mass; Column 
3, the virial radius; Column 4, the radial extent of the 
gas disk, defined as the largest radius at which we find 
gas particles in the disk plane. Columns 5 and 6 contain 
the mass of each gas and DM particle in the highest res- 
olution region, and Columns 7 and 8 are the softening 
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TABLE 1 
Properties of simulations 



Name 


M vir 


^vir 


^disk 




^DM 


e g as 




Qo ho 






(M ) 


(kpc) 


(kpc) 


(M ) 


(M ) 


(kpc) 


(kpc) 






KGCD 


8.8 x 10 11 


240 


10 


9.2 x 10 5 


6.2 x 10 6 


0.57 


1.1 


0.3 0.7 


0.039 


AGCD 


9.3 x 10 11 


270 


21 


3.3 x 10 6 


1.9 x 10 7 


0.87 


1.5 


0.3 0.65 


0.045 
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Fig. 1. — Hi gas column density maps of KGCD (left set of 4 
panels) and AGCD (right set) simulations. 




lengths in that region. The cosmological parameters are 
presented in Columns 9-11. Note that the spatial resolu- 
tion for the gas is determined by the smoothing length of 
the smoothed particle hydrodynamics scheme. The min- 
imum smoothing length is set to be half of the softenin g 
length of the gas particles fsee lKawata fc Gibsonl l2003). 
The smoothing length depends on the density, and the 
average smoothing length in the simulated HVCs we fo- 
cus on here is ~ 5 kpc. 

Both galaxies are similar in size and mass to the MW, 
and have clear gas and stellar disk components. Figure E] 
shows edge-on and face-on views of the projected gas 
density of each galaxy at the final timestep. We use the 
simulation output at z = 0.1 for KGCD, as contamina- 
tion from low-resolution particles in the simulated galaxy 
start to become significant at this redshift. We use the 
output at z = for AGCD. 

In order to compare the simulations with the HVC ob- 
servations in the MW, we set the position of the "ob- 
server" to an arbitrary position on the disk plane of 
the simulated galaxies, with galactocentric distance of 
8.5 kpc, and "observe" the Hi column density and ve- 
locity of the gas particles from that position. Figure [2| 
demonstrates the H I column density of HVCs using all- 
sky Aitoff projections. Here, we define HVCs as consist- 
ing of gas particles with absolute lines-of-sight velocities, 
f LSRo deviating from the Local Standard of Rest (LSR) 
by more than 100 km s _1 — Figure[2] shows the H I column 
density of these gas particles with |^lsr| > 100 km s -1 . 
We set the rotation velocit y of the LSR to 220km s _1 , 
similar to that of the MW lLockman et all (|2002ft . We 
can confirm that both simulated galaxies have gas disks 
rotating at ~ 220 kms -1 . In this paper, results are based 
only on those particles within two virial radii (r v i r ; see 
Table 0). We have also confirmed that the results are not 
sensitive to the cutoff radius chosen for column densities 




Fig. 2. — Hi column density maps of KGCD (top) and AGCD 
(bottom) simulations, in a all-sky Aitoff projection in using galac- 
tic coordinates. Particles with I^lsrI < 10 kms -1 are excluded. 
Compare these maps with fig. 2 in IWakkeri |1991). 



7V(HI) > 10 1T cm -2 . We only display results for one 
chosen observer, however we have confirmed the general- 
ity of these results, with the sky coverage fraction typi- 
cally changing by no more than 20% for a given column 
density, as we change the observer's position and/or we 
analyze other outputs of the simulation near the final 
redshift. 

Our chemodynamical simulation follows the hydrogen 
and other elemental abundances for each gas particle, 
but does not calculate the ionization fraction of each 
species. Instead, the Hi mass for each gas particle 
is calculated assuming collisional ionization equilibrium 
(CIE). The CIE neutral hydrogen fra ction is estimated 
using CLOUDY94 (Ferlan det alJ ll998). We multiply the 
fraction by the hydrogen abundance, and obtain a Hi 
mass fraction for each particle as a function of its tem- 
perature. We ignore any effect from the background ra- 
diation field. 

3. RESULTS AND DISCUSSION 

Figure [3 demonstrates that both galaxies have a signif- 
icant number of HVCs, with column d ensities compara- 
ble to those observed bv lWakkerl (|1991f) . KGCD displays 
several large linear HVCs at galactic longitudes / ~ 60° 
and / ~ 270°. These components correspond to the outer 
ring structure seen at a galactocentric radius of ~ 30 kpc 
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Fig. 3. — Hi sky coverage fraction vs. limiting column 
density. The dashed line denotes the sky coverage from the 
ILockman et alJ 120021) HVC survey (excluding velocities ^lsr i n 
the range ±100 km s -1 , and with poor coverage of latitudes \b\ < 
20°) to the 4cr completeness limit. The sky coverage obtained from 
the KGCD simulation is shown as a solid line, and the dotted line 
denotes the AGCD simulation. 



in Figure^ We name this structure the "polar gas ring" , 
and discuss it later. 

To compare the HVC population of our simulations 
with the MW HVCs quantitatively, Figure shows the 
fraction of sky covered by HVCs as a function of limiting 
column density for both simula tions and for the observa- 
tions in ILockman et all ((2002). In this plot, we exclude 
the area with low galactic latitude \b\ < 20°, to avoid 
contaminatio n by the disk com ponent (the sample of 
sightlines in ILockman et all200 2 was limited in a similar 
fashion). As is obvious from Figure [2j KGCD has more 
high column density HVCs than AGCD, and almost all 
of the sky is covered down to 10 16 cm -2 . At a fixed col- 
umn density, the sky coverage in the simulations bracket 
the observed sky coverage in the MW. Note that we ig- 
nore any effects of background radiation. It is expected 
that such a field would decrease the population of HVCs 
with Galactocentric distanc e less than 10 kpc or N (HI) < 
10 19 cm- 2 (|MalonevllT993 iBland-H awtho rn fc Malonevl 
[T9971 [l99l . Thus, since as mentioned below, the dis- 
tances of the simulated HVCs are greater than 10 kpc, 
our estimated coverage fractions should be interpreted as 
an upper limit for 7V(HI) < 10 19 cm -2 , and a lower limit 
for 7V(HI) > 10 19 cm -2 , where the highest column den- 
sity HVCs may not be fully resolved. With these caveats, 
KGCD appears to have a sufficient population of HVCs 
to explain the observed population of HVCs within the 
MW. We conclude that current cosmological simulations 
can produce MW-size disk galaxies with similar popu- 
lations of HVCs to those in the MW. The differences 
between KGCD and AGCD may demonstrate real differ- 
ences in the populations of HVCs among disk galaxies. 
However, to understand the causes of such differences, 
we need a larger sample of high-resolution simulated disk 
galaxies, the subject of a future study. 

We find that velocities of the simulated HVCs seen 
in Figure El are d i stribu ted similarly to those observed 
bv iPutman et all ((20021) . Overall, the clouds between 
I = 0° and 180° have a negative velocity, while the 
clouds between I = 180° and 360° have a positive ve- 
locity. This is natural, since the LSR is moving to- 
wards I = 90°. In KGCD, we find that there is a rel- 
atively large HVC complex whose velocity is very high, 
—450 < i>lsr ^ — 300kms _1 . These very high velocity 
clouds (VHVCs) are located between -45° < b < 15° 
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Fig. 4. — Hi first moment (column density weighted mean veloc- 
ity in the LSR) map of KGCD (top) and AGCD (bottom) simula- 
tions. Particles with |i>ls r| < 100 km s" 1 are excluded. Compare 
these maps with fig. 16 in Wakker et al. (2003). 



and 60° < I < 120°, in the top panel of Figure H 
The MW also has such VHVCs in the Anti-Cente r com- 
plex ( Hu lsboschllT973 iHul sbosch & Wakkejl ll98fih . The 
galactocentric distance to the VHVC in KGCD is ~ 10- 
25 kpc, and we find that the cloud is a gas clump which 
has recently fallen into the galaxy. We convert the 
Hi mass for e ach particle into a Hi flux us ing Mhi = 
0.235^ pc 5 t ot (|Wakker fc van Woerdenl 19911) . where the 
Hi mass Mhi is in M , total Hi flux in Jy kms -1 is 
5totj an d distance from observer to the particle in kpc 
is £>k P c- We find the Hi mass of the VHVC complex 
is 1.2 x 1O 8 M , and the total Hi flux (for the chosen 
observer) is 2.5 x 10 6 Jy kms -1 . Since it is infalling on 
a retrograde orbit, its velocity relative to the LSR be- 
comes very large depending on its location relative to 
the observer. Therefore, the observed VHVCs may be 
explained by such infalling gas clumps within the Galaxy. 
It is also worth noting that we do not find any associated 
stellar or DM components in the simulated VHVCs. 

In the simulations, we are able to measure the distance 
of HVCs — data which are not yet generally available for 
the real MW. In Figure [5J we show a flux- weighted his- 
togram of the galactocentric distances, tmw, of the high 
velocity gas particles with \b\ > 20°. In both simula- 
tions, less than 1% of the Hi flux visible in Figure 
originates from HVCs with tmw > 150 kpc. This is 
consistent wi t h the aforementioned limits established by 
iPisano et all (|2004l) . 

It is also clear that most of the emission in the KGCD 
simulation results from the polar gas ring, whose radius 
is ~ 30 kpc. The mass of this ring (including low veloc- 
ity gas) is 3.6 x 1O 9 M , of which 2.3 x 10 9 M is Hi. 
Although there is no prominent polar ring in AGCD, at 
larger distances, the flux distribution is found to follow 
that of the KGCD simulation. 

Figure demonstrates that the polar gas ring forms a 
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Fig. 5. — Hi flux weighted histogram with logarithmically spaced 
bins in distance (de-emphasing fluxes from smaller galactocen- 
tric distances) of both AGCD and KGCD simulations, excluding 
I^lsrI < 100 krns -1 and \b\ < 20°. 



linear high velocity structure found in all quadrants of 
the sky. Large HVC components, such as Complex C 
and the Magellanic Strea m, are well kn own in our MW, 
and several authors (e.g. iHaudl [l988) argue that the 
MW is surrounded by a polar gas ring. Neglecting the 
Magellanic Stream, which ha s likely originated from the 
infalling Magellanic c louds (|Ga^iner^^^offuclril 1996; 
lYoshizawa fc Noguchi \200A IConnors et all 120050 . the 
largest HVC structure is Complex C. We measure the 
mass of the HVC ring in one quadrant as ~ 4 x 10 8 M . 
If we place Complex C at a galactocentric distance of 
30 kpc (only a lower limit has been obtained for its galac - 
tocentric distance of 8.8 kpc; Ivan Woerden et aLl ll999). 
its mass would be ~ 1.5 x 1O T M . Nevertheless, this 
indicates that there is enough H I gas to create such large 
observed HVCs in the simulated galaxy. In KGCD, the 
polar gas ring is a relatively recently formed structure 
that begins forming at redshift z ~ 0.2, and prior to this 
time, the associated gas particles are found flowing in- 
ward along filamentary structures. Thus, this simulated 
ring structure demonstrates that current ACDM numer- 
ical simulations can explain the existence of such large 
HVCs like Complex C, as recently accreted gas which ro- 
tates in a near-polar orbit. On the other hand, AGCD 
does not show such a prominent large structure. This 
may indicate that such HVC structures are not common 
in all disk galaxies. Although this study focuses on only 
two simulated galaxies due to the limit of our computa- 
tional resources, more samples of high-resolution simu- 
lated galaxies would elucidate how common such large 
structures are, and what kind of evolution history is re- 
quired to make such HVCs. 

We also analyze the metallicity of the simulated HVCs. 
In KGCD, we find that the prominent HVCs have Hi 
flux weighted metallicities of —4 < \og(Z/ Z ) < —2 
with a flux weighted mean of \og(Z / Z ) ~ —2.4. This is 
much lower than the metallicities of the observed HVC s 
in the MW, -1 < log(Z/Z0) < (iWakkerJ Ejool . 
The polar gas ring seen in KGCD has a metallicity of 
—3 < log(Z/Zo) < 0.5 with a mean of log(Z/Z ) ~ 
— 1.7, which is lower than the obs erved metallicity of 
Complex C, logfZ/ Zp) - -1 (e.g. iGihson eTaTl f2001: 
Rich ter et alJ 120011) . Thus, our numerical simulations 
seem to underestimate the metallicity of the later in- 
falling gas clouds. This is likely because we adopt a weak 
supernova feedback model in our simulations. If we use a 
model with strong feedback, more enriched gas is blown 
out from the system at high redshift, which can enrich 



the inter-galactic medium which then falls into galaxies 
at a later epoch. 

The majority of the simulated HVCs, including the 
polar gas ring, do not have any obviously associated 
stellar or DM comp onents, which is consistent with 
current observations (Simon & Blitz 200^1 iSieerel et alJ 
2005). However, a few compact HVCs are found to be 
associated with stellar components. It would be inter- 
esting to estimate how bright they are, and if they are 
detectable within the current observational limits. Un- 
fortunately, the resolution of the current simulations are 
too poor to estimate the luminosity, and it is also likely 
that our simulations produce too many stars due to our 
assumed minimal effect of supernova feedback. 

The clouds in the simulation may be destroyed by ef- 
fects that our simulations are not able to reproduce accu- 
rately. The resolution and nature of the SPH simulations 
make it difficult to resolve shocks between the HVCs 
and the MW halo; however simple analytic estimates 
show that there will not be strong s hocking of the HVCs 
due to the low density of the halo. Ilndebetouw fc Shull 
(2004), using more detailed simulations, argue that the 
Mach number of the HVCs are only M s ~ 1.2-1.5, 
marginally sufficient to form s hocks, and heat i ng the 
leading ~ 0.1 pc of the HVC. Qnilis Moorel GflS3) 
show that the lifetime of such HVCs in the presence of 
shocks is rsj 1 Gyr, long e nough for our HVCs to survive. 
Mail er fc Bullockl (|2004l) discuss various physical pro- 
cesses, including conduction, evaporation, ram-pressure 
drag, Jeans instabilities, and Kelvin-Helmholtz instabili- 
ties that limit the mass ranges of stable clouds. The most 
massive HVC apparent in our simulations is found to be 
1.5 x 1O 8 M , and the mass resolution and hence small- 
est resolvable HVC in our simulations are ~ 1O 6 M ; 
clouds of both extremes are cle arly within the stable mas s 
ranges summarized in fig. 6 of lMaller fc Bullockl ((2004). 
Therefore, the simulated clouds are also expected to be 
stable. 

We report that HVCs seem to be a natural occurrence 
in a ACDM Universe. We emphasize that the galaxies 
that result from our simulations were not created specif- 
ically to reproduce the MW exactly — they were selected 
for resimulation at higher resolution on the basis of being 
disk-like L* galaxies. However, we have serendipitously 
discovered that simulated galaxies that are similar in size 
to the MW naturally contain H I gas in the vicinity of the 
disk that are similar to the anomalous velocity features 
seen in the MW. 
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